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Abstract. 

Surfaces in i-Al e8 Pd 2 3Mn 9 as observed with STM and LEED experiments 
show atomic terraces in a Fibonacci spacing. We analyze them in a bulk 
tiling model due to Elser which incorporates many experimental data. The 
model has dodecahedral Bergman clusters within an icosahedral tiling T*( 2F ) 
and is projected from the 6D face-centered hypercubic lattice. We derive 
the occurrence and Fibonacci spacing of atomic planes perpendicular to any 
5fold axis, compute the variation of planar atomic densities, and determine 
the (auto-) correlation functions. Upon interpreting the planes as terraces at 
the surface we find quantitative agreement with the STM experiments. 

1 Introduction. 



The bulk structure of the icosahedral phases i-AlPdMn, i-AlFeCu and their 
modelling in terms of a 6D description has been an active research field for 
more than one decade. From the many papers on this topic we mention only 
a few. Many more references are quoted in these publications. Katz and 
Gratias [|7| derive from previous work for i-AlFeCu a quasiperiodic network 
of atomic positions. It is generated by three basic atomic windows related 
to the 6D hypercubic F-lattice. They examine carefully the interatomic 
distances. De Boissieu et al. (2| determine for i-AlPdMn from X-ray and 
neutron data in detail the decomposition of the atomic surfaces. All these 
models use the 6D embedding, the parallel and the perpendicular space. 
Elser m generalizes and unifies both these models and interprets them in 
terms of clusters occupying the odd and even vertices of the icosahedral tiling 
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related to the 6D hypercubic P-lattice: The odd vertices form the centers 
of Bergman clusters, which then around the even vertices build up Mackay 
clusters. Additional atomic positions are related to this basic structure. The 
Elser model was actually created for the study of random tilings, but by 
construction admits a perfect tiling structure which then incorporates the 
main experimental data which have led to the models by Katz and Gratias 
and by de Boissieu et al. As shown in [pTH , the Elser model can be taken as a 
network of atomic positions in a tiling model, denoted by T*( 2F ^ and related 
to the observed hypercubic F-lattice and -module. This tiling model will be 
used in what follows. Its composite atomic surfaces are closely related to 
those of the Katz-Gratias model 0. 

The surface structure of i-Al§%Pd22,Mn$ perpendicular to 5fold axes has been 
explored by various groups. Schaub et al. |15| applied scanning tunneling 
microscopy (STM) and low-energy electron diffraction (LEED) to obtain 
atomic scale information of a sputtered and annealed surface. They observe 
a sequence of 11 atomically flat terraces. Two spacings of (4.22 ± 0.26) 
and (6.78 ± 0.24) A form a Fibonacci string of the type LLSLLSLSLL. 
Pentagonal holes of a single fixed orientation appear within these terraces. 
Gierer et al f| ||] by dynamical LEED studies for a similarly prepared surface 
confirm the quasicrystalline structure. To interpret their data they perform 
dynamical diffraction calculations for assumed terminations of a model patch 
from the bulk model of de Boissieu et al. (2|]. They find optimal agreement 
for Al-rich terminations of high atomic density. A study by Ebert et al. of 
in-situ cleaved surfaces revealed terraces only after annealing of the initially 
rather rough surface. 

For the theory of quasicrystals, the experiments raise the question what 
quasiperiodic repetition pattern and what structure variation within planes 
can be rigorously obtained from a bulk model of i-A/ggPG^Mng. In the 
present paper we address these questions in the description by the tiling 
model. We develop a quasiperiodic analysis similar to the one of crystal sur- 
faces in terms of particular netplanes. To obtain exact results we cannot rely 
on features seen in a model patch. Instead we make full use of the technique 
of windows or coding for quasiperiodic structures. The general principles of 
the window technique have been described in several monographs on qua- 
sicrystals. We illustrate them on the well-known one-dimensional Fibonacci 
paradigm. We then apply the unique lifting and projection method between 
the physical (parallel) and the window (perpendicular) space, called the star- 
map by Moody |14]| , to the icosahedral F-module, to the tiling, and to its 
decoration. It turns out that our main results can be expressed in terms of 
the one-dimensional Fibonacci system. 
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We now survey the model input and the content of the following sections. 
Our bare tiling model has the following data: We adopt the T*( 2F ^ tiling 
model projected from the face-centered hypercubic lattice (2F) ~ D e in E 6 . 
Upon scaling by a factor 2, the lattice (2F) comprises the even vertex points 
(even index sum) of the full hypercubic lattice P whose projection was given 



in 0. For a full description of the tiling and its projection we refer to pOfl . Its 
vertex points are projected lattice points. We shall use two units of length: 
© is the length along 5fold axes of the six basis vectors e i? i = 1, . . . , 6 of 
the hypercubic lattice, projected to the two invariant icosahedral subspaces 
E\\ or E± respectively. Along projected 2fold axes we choose the standard 
length (2) = ^7=^j©- To convert to atomic distances in i-AlPdMn we adopt 
-scaled short edge length 

©^4.56 A. 



from ||, |TTJ the r-scaled short edge length of the tiling 

2r 



The window of the vertex points for the tiling is in E± the triacontahedron 



||, [T0|] shown in Fig. 5. The tiling is decorated according to Elser || with 
dodecahedral Bergman clusters [Tl|| . The midpoints of these Bergman clus- 
ters are placed on the projected odd vertex points of the hypercubic lattice. 
Their edge length is r -1 © = 2.96 A , their height along a 5fold direction is 
© = 6-60 A . For all other atomic positions, most of which do not enter 
the present analysis, we refer to [|, 

In section 2 we develop the window technique for the bulk tiling and its planes 
perpendicular to 5fold axes. We start in 2.1 with the Fibonacci tiling and 
explain the technique of windows. We briefly describe the icosahedral tiling 
q-*(2F) f Qr F_ph ase) section 2.2, and planes of vertex points perpendicular 
to a 5fold axis in a 3D space E\\, section 2.3, and give their windows in E±. 
We shift between these planes along Fibonacci lines, section 2.4, and show 
that most of the vertex points belong to a system of shifted planes, section 
2.5. In section 2.6 we interpret the terraces found in the STM experiment as 
terminations of the bulk model. From the tiling model we prove the existence 
of a full Fibonacci sequence of planes and of a spacing as found in the STM 
experiment, and we predict variations of the density of vertex points along 
the sequence, with bounds from the observed Fibonacci string. 

In section 3 we use the decoration of the tiling to infer more structure infor- 
mation within the planes from other atomic positions of the tiling model. In 
particular we look for pentagonal structures as seen in the STM experiments 
|To| . We consider the dodecahedral Bergman clusters of the Elser model |^| 
on the tiling. The dodecahedra have two pentagonal vertex sets of the same 
orientation perpendicular to a 5fold axis. The corresponding cutting planes 



3 



are transformed in subsection 3.1 by lifting and projection to the perpendic- 
ular space. Their window description with respect to the triacontahedron is 
derived. The correlation with vertex points of the tiling gives rise to three 
alternative models for the structure within planes. The predicted density of 
vertex points and pentagons is derived in exact form in section 3.2. In sub- 
section 3.3 we compute in closed form the Patterson function within planes 
for vertex points and pentagon centers. 

The bulk structure of the tiling model, analyzed here up to the level of 
Bergman clusters, displays for the planes a repetition and structure pattern in 
line with the terrace structure found in the experiments [n| which stimulated 
the present analysis. A complementary approach to the terrace structure, 
based on generating a model patch, is given in [|1^| and confirms the present 
analysis. 

2 Tilings and windows. 

2.1 Fibonacci lines, their windows, and search for the 
string LLSLLSLSLL. 

We recall the well-known projection and window technique for the Fibonacci 
tiling T. We shall emphasize the window technique since it will be needed 
when we apply in sections 2.4-2.6 Fibonacci lines to the icosahedral tiling. 
Let A be the square lattice in 2D whose edge length we adjust for convenience 
to a/t + 2, t = |(1 + In a lattice basis the points of A are 

x = n\e\ + n 2 e 2 . (2) 

In a system of coordinates (x\\,x±) rotated by : c = cos(0) = ~^| , s = 
sin(0) = -^g wrt. the natural basis, the basis vectors are 

d = (c, -s)Vt + 2, e 2 = (s, c)Vr + 2, (3) 

and the coordinates of the lattice points become 

(x\\,x±) = (nidn + n 2 e 2 \\,n 1 e 1± + n 2 e 2± ) (4) 
= (x||(ni,n 2 ),x_L( , n-i,^2)) = {n\T + n 2) —n\ + n 2 r). 

The projections (x\\,x±) form two r-modules on orthogonal lines E\\,E± re- 
spectively. There is a unique map x\\(ni,n 2 ) ^ x±(rii,n 2 ) between these 
modules, corresponding to the star map of Moody HTH , and there is a unique 
lifting of x\\ {ni,n 2 ) or x±(jix, n 2 ) into a point of A. The projections x\\ (A), x±(A) 
cover Eh , E± dense and uniformly. 
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Upon choosing in E± the window f± := (—1, r], whose length \w±_\ = r + 1 is 
the projection of a unit square to E±, the vertex set of the Fibonacci tiling 
T in E\\ is defined as 

v(T) = {xi\(n 1 ,n 2 )\x ± (n 1 ,n 2 ) G f±}. (5) 

The endpoints of the window are restricted in order to avoid ambiguities. 
When T is lifted back into A e E 2 , it forms the vertex set of a continuous 
staircase formed by edge lines as shown in Fig.l. 
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Fig.l The Fibonacci tiling is the projection of a staircase, formed by edge 
lines in a square lattice, to a line E\\ of slope r -1 . The successive vertices x\\ of 
the staircase may be enumerated by the single integer N = ni+n 2 . Projected 
to x\\, the edge lines form the Fibonacci tiling with two tiles S,L of length 1, r 
respectively. The projections x±(N) of the vertices to the orthogonal space 
E± fall into a window f± = (— 1, r] of length r + 1. 

The projections of the steps to E\\ form the familiar Fibonacci tiling with 
two tiles S,T of length l,r respectively. With x\\ increasing, adjacent tiles 
form the vertex configurations LS, LL, or LS respectively. The windows in 
E± for these vertex configurations can be shown to form subwindows of f± 
given by 

("1,0], (6) 
(0,r-l], 
(r-l,r]. 

We now wish to compare and analyze Fibonacci tilings with different starting 
points. Because of the uniform dense covering, we may choose in E± an 



fLS 
J± 
f LL 
J± 
i-SL 
J± 
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arbitrary initial point c± G w± and associate to it an initial point of a tiling 
T(cjJ. We label the initial vertex by (0,0) =>- and the successive vertices 
of T(c_i_) by the single integer N = n\ + rii- From the window condition we 
can generate x\\(N), x±(N) step by step according to 



x±(N + l) 
xn(N + l) 



x±(N) - 1 ^ (x ± (N) - 1) e f ± 
x ± (N) + T^(x ± (N) + T)ef ± 

X\\(N)+T, 

x«(N) + l. 



(7) 



The steps in E\\ propagate the tiling by a new tile L or S respectively. For 
later purposes, like the determination of densities of points in section 3.2, we 
emphasize the propagation as a function of N in terms of the perpendicular 
coordinate in the window. 




0123456789 101112 131415 16171819 20 2122 23 24 



Fig. 2 Four Fibonacci lines starting at an LL vertex are coded by four initial 
points in a vertical subwindow scaled by r~ 3 . The vertical coordinate is 
y_i_(N) of eq.(8). For steps numbered from to 24, the lines connect the 
images in the window for these four points. Each step produces in E\\ a long 
or short interval of the corresponding Fibonacci line. 

We adjust the perpendicular coordinate to the midpoint of the window and 
scale it by a factor r to obtain the new variable 

y±(N) := rx ± (N) - \ (8) 
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whose window w± = (— |t 3 , ^t 3 now has the length \w±\ = r 3 with the 
central subwindow for LL vertex configurations of length |w^ L | = 1. It can 
be shown that the subset of LL vertices of the original tiling form another 
Fibonacci vertex set scaled by a factor r 3 . The function y±(N) is plotted in 
Fig. 2 for four initial values from the subwindow w^ 1 . This Figure illustrates 
the variety of sequences as a function of the initial value. Successive values 
are connected by straight lines. The reason for starting at an LL subwindow 
will become apparent when we go to the icosahedral tiling in section 2.4. 
With the window technique we search for the finite string LLSLLSLSLL 



found in the terrace spacing of the experiment p5 |. For a Fibonacci line 
coded by the initial point y±(0) = — \, this string occurs at the points N= 9 
. . . 19, compare Figs. 2,6. For other initial points, the string would occur at 
some other step. We infer all possible occurrences of the string as conditions 
with respect to the window: The string will be stable under vertical shifts 
Az/_l(0) of the initial point as long as its highest value y_i_(17) and its lowest 
value 2/j_ (14) do not pass the limits ±|r 3 respectively of the window w±. 
These window conditions are independent of the initial point. Clearly the 
appearance of the string puts narrow bounds on the corresponding values of 
y±, compare section 2.6. 



2.2 Icosahedral tilings. 

The construction of 3D tilings follows the paradigm given by the Fibonacci 
line. The projections are now determined by requiring non-crystallographic 
and in particular icosahedral point symmetry after projection. It is well- 
known that an icosahedral tiling T p with two rhombohedral tiles arises by 
icosahedral projection to 3D from the primitive hypercubic P-lattice and 
module in 6D 0. 

In two orthogonal 3D spaces E\\, E± we find the six 5fold, ten 3fold and 
fifteen 2fold axes associated with the icosahedral group. The six primitive 
basis vector ei, . . . ee of the hypercubic lattice upon projection point along 
5fold axes. Their length we denote by (5), and their directions we choose 
as follows [|i~0| : In E\\ we take cos(ei||, e*|| ) = l/y/r + 2, i = 2 ... 6 and for 



i 



2. . .5 pass from e^i to e i+ i\\ by a rotation around en with angle 27r/5. 
In E± we take cos(eu_, e^jj = —l/y/r + 2, i = 2 . . . 6 and for % — 2 ... 5 pass 
from ei± to e i+ i± by a rotation around ei± with angle 4n/5. All vectors along 
2fold axes arise from projections of (e^ ± e 3 -), i ^ j. Their shorter length we 
denote by © = ^=©. 

The icosahedral quasicrystals i-AlFeCu and i-AlPdMn from their diffraction 
pattern are indexed by the hypercubic face-centered or F-lattice and -module 
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rather than the primitive P-module. The hypercubic F-lattice, scaled by a 
factor 2 and denoted here as (2F), may be viewed as the subset of even lattice 
points (even index sum) from the full hypercubic lattice in 6D. 
Turning attention to this lattice and module, we briefly summarize the con- 
struction of the icosahedral tiling T*( 2F ) associated with the F-lattice and 



given in [10] . In both 3D spaces we have 6D modules whose bases can be 
formed for example from three short and three long vectors along three se- 
lected 2fold axes. By the unique lifting and projection, there is a one-to-one 
map q\\ <=>■ q±, the star map of [fy[], between points q\\ and q± of the two 
modules in E\\ and E±. For simplicity we suppress the basis and the six 
integers in q\\,q± which generalize eq.(4) and underlie this map. The vertex 
points of the tiling T*" > generalization of eq.(5), given by 

v(T*( 2F >) = {q\\\q± G triacontahedron}, (9) 

i.e. the projections q\\ G E\\ of those lattice points whose projections q± G E± 
fall into the triacontahedral window, compare Fig. 5. The projections {q±} 
fill the triacontahedron dense and uniformly. The triacontahedron is the 
icosahedral projection to E± of the Voronoi or Wigner-Seitz cell of the F- 
lattice in 6D. The tiling T*^ 2 ^ has six tetrahedral tiles. In its present simple 
form we need only two tetrahedra with 3fold symmetry axis. The vertices of 
these two tetrahedra coincide with four even vertices of the two rhombohedra 
associated with the tiling T p . The simple form of the tiling T*( 2F > is fully 
described by putting atoms into positions on the full rhombohedral tiles but 
allowing for the distinction of even and odd vertices, as is done in the Elser 



model U|, pT| . We shall need only the even and odd vertex points. 
The relation of the tilings T*( 2F ) and T p may be summarized as follows: The 
triacontahedral windows for the vertex sets coincide. The modules differ from 
one another: The (2F) module is the even submodule of the P module. By 
expanding each of the two 3fold symmetric tetrahedra back into the corre- 
sponding rhombohedron and dropping the distinction between even and odd 
vertices we can locally derive T p from T*^ 2F \ 



2.3 Planes perpendicular to 5fold axes. 

We turn to planes in the tiling T*( 2F \ Fix in E\\ a 5fold axis parallel to e\\\ 
as in Fig. 4 and consider vertex points q\\ in a plane perpendicular to it. Next 
we pass to E±, consider the corresponding 5fold axis parallel to e%± and the 
images q± of the vertex points q\\ from the plane under the one-to-one map. 
It turns out that these images q± G Ej_ lie again in a plane perpendicular 
to the 5fold axis. In addition they must be points from the triacontahedron. 
Hence we get the result: 
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Fig. 3 Kepler's decomposition of the triacontahedron into a central decagonal 
prism Xx and two shells. 

Prop 1: The window for vertex points q\\ from the tiling T*( 2F ) in a fixed 
plane perpendicular to a 5fold axis is in E± the intersection of a plane per- 
pendicular to the corresponding 5fold axis with the triacontahedron. 
The triacontahedron is shown in Fig. 5 in a view perpendicular to a 5fold 
axis. The distance from the center to a 5fold vertex is r (5), where © is the 
standard length along a 5fold axis. The triacontahedron with respect to this 
5fold axis has a central decagonal prism of thickness ^-r (5). J. Kepler in 
1619 || not only introduced the triacontahedron, but also visualized these 
decagonal prisms and denoted them by the letters Xx, see Fig. 3. 
The central decagonal prism of the triacontahedron when seen as a subwin- 
dow for part of the tiling 7"*( 2F ) has a particular significance, as is shown in 
fli~3"|| : Any planar decagonal intersection of the triacontahedron in E± deter- 
mines in Eh an infinite planar tiling TTT by two golden triangles, compare 
With respect to the full 3D tiling T*( 2F \ this planar subtiling is formed 
by faces of tetrahedral tiles. 

The planar intersections of the triacontahedron outside the decagonal prism 
are windows for planes of vertex points which in general do not form a planar 
tiling. From the uniform covering of the window it follows that the density 
of vertex points in all planes is proportional to the area of its window, i.e. 
of the corresponding intersection of the triacontahedron. We shall compute 
this density in subsection 3.2. 
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2.4 Fibonacci shifts between parallel planes. 



The planar TTT subtiling has the property that through any vertex point 
there passes at least one infinite Fibonacci line. In terms of its decagonal 
subwindow this results from the geometric property that any interior point 
belongs to at least one subwindow for an infinite Fibonacci line. In E\\ this 
Fibonacci line points along a 2fold axis associated with two vectors whose 
length scales by r. All subwindows for a fixed Fibonacci line are sections 
of length t 2 (2) on parallel lines perpendicular to and bounded by opposite 
rectangular faces of a decagonal prism. 



2' 




Fig. 4 Two typical perpendicular 2fold axes 2, 2' can be chosen within a plane 
with a 5fold axis 5. A plane perpendicular to this axis can be shifted from 
the origin by vectors along these 2fold axes. 

As the initial plane we shall choose a reference plane whose vertex points 
form a triangle TTT pattern. All these planes have the same highest density 
of vertex points, see subsection 3.2. To shift between planes perpendicular 
to a fixed 5fold axis, we shall use vectors along 2fold axes outside this plane. 
From the orbit in E\\ of 2fold axis with respect to the 5fold one we pick two 
perpendicular 2fold coplanar axes 2, 2' which form with the 5axis the angles 
arccos(-^=) = 58.3 degrees, arccos(-^=|) = 31.7 degrees, see Fig. 4. In the 
notation of [|l(|, we choose the axis 5 along en, the axis 2 along the short 



and long vectors — (e 2 + e 3 )|i, (ei + e 5 )||, and the axis 2' along the short and 
long vectors (ei + es) j| , — (e^ + eo)\\ respectively. 

The vectors along the 2fold axis 2 from eq.(l) have the r-scaled model length 
r (2), r 2 (2) . By multiplication with the cosine of the corresponding angle we 
get the parallel spacings of planes perpendicular to the 5fold axis 5. These 
spacings become ©) and {^r^ ©) respectively. By comparison, the 
two vectors along the 2fold axis 2' yield along the 5fold axis spacings scaled by 
a factor r. Therefore the latter vectors will not generate additional parallel 
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planes. We obtain the short and long spacings 4.08 and 6.60 A respectively, 



fully in line with the terrace spacing observed in ||15|| and quoted in section 1. 
So we have identified in the bulk tiling model the shift vectors which generate 
the terrace structure. 




2' 



Fig. 5 Three decagonal prisms of the triacontahedron in E±: The first one 
has its 5fold axis 5 in the vertical direction. The second and third prism have 
opposite rectangular faces perpendicular to two 2fold axes 2,2' coplanar with 
the axis 5. 

To assure that, starting from a fixed plane, we generate by shifts an infinite 
system of parallel planes, we turn to E±. In E± the three axes 2, 2', 5 remain 
coplanar but the directions and angles of 2 and 2' are interchanged, see Fig. 5. 
We get infinite Fibonacci lines along the axes 2,2' if the vectors along these 
axes can be associated to decagonal prisms of the triacontahedron. Two 
decagonal prisms with this property and associated with 2,2' respectively are 
shown in Fig. 5. 

The initial reference plane was chosen with the triangle tiling and hence has 
as its window a decagonal section through the triacontahedron perpendicular 
to the 5fold axis. An infinite system of parallel planes will arise, provided that 
we select a starting point which also belongs to a Fibonacci window associated 
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with vectors say along the axis 2. The projection of the full Fibonacci window 
along the 5fold axis in E± by multiplication with cos(2, 5) = ^l +2 becomes 

T^Tj©- The projection of its central subwindow for LL vertices equals the 
thickness -^-^ © of the decagonal prism. We conclude that, among the parallel 
planes shifted along the infinite Fibonacci line in the direction 2, the reference 
plane and in fact any dense plane occurs at the LL vertices. 

We summarize the information obtained so far on parallel planes of vertex 
points q\\ in the tiling T*^ 2F \ generated by Fibonacci lines along 2fold axes: 
Starting from a dense reference plane, we generate an infinite set of parallel 
planes. They follow a Fibonacci spacing with perpendicular short and long 
distances © = 4.08 A and ^ © = 6.60 A, fully in line with the STM 
observations of terraces. The dense planes occur at all LL vertices of the 
generating Fibonacci line. Other parallel planes in the set will have a lower 
density of vertex points. The string LLSLLSLSLL analyzed in section 2.1 
can now be converted into a sequence of parallel planes of varying density. 
To complete the analysis of parallel planes, we must find out what fraction 
of all vertex points q\\ E T*^ 2 ^ is reached within this infinite sequence of 
planes. 

2.5 Parallel planes are connected by Fibonacci lines. 

We wish to show that indeed we can reach from a fixed dense reference plane 
most vertex points q\\ by shifts along Fibonacci lines. For this purpose we 
consider only those vertex points q\\ which lie on at least one infinite Fibonacci 
line. From the window side we know that this is the case if q± is a point 
from any decagonal prism. There are six such prisms, and this motivates the 
definition of a new window: 

Def 2: Decagonal prism Xx approximation: We analyse only those points q± 
of the triacontahedral window which belong to at least one decagonal prism, 
that is to the union prism j. We omit in this approximation the points q± 
from small parts of the triacontahedron close to its 5fold vertices, compare 
Fig.5. 

The vertex points q\\, with q± belonging to this new window, have the fol- 
lowing properties: 

Prop 3: Consider vertex points q\\ in a plane parallel to a fixed dense TTT 
infinite reference plane. Among them there is a point on an infinite Fibonacci 
line which intersects (as a continuous line) the reference plane. 
Proof. Through any point coded in the decagonal prism approximation there 
runs at least one infinite Fibonacci line. If it intersects the reference plane 
we are through. If it runs parallel to the reference plane, we can (proof omit- 
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ted) in at most two parallel steps pass to another point with an intersecting 
Fibonacci line. 

Prop 4: If an infinite Fibonacci line intersects as a continous line a TTT 
plane, it hits this plane in a vertex point. 

Proof. The points of the infinite reference plane form the vertices of the planar 
TTT subtiling by faces, the points of the non-parallel infinite Fibonacci line 
form the vertices of a linear subtiling by edges of the 3D tiling T*^ 2F \ Both 
subtilings are parts, hence their intersection is a vertex of the full tiling. 

The two propositions allow us to code in E± planes of vertex points, parallel 
in £j| to a first dense reference plane perpendicular to 5, by their intersections 
with Fibonacci lines along the axis 2 of Fig. 4: 

Prop 5: Any vertex plane perpendicular to a 5fold axis has at least one point 
connected to the reference plane by an infinite Fibonacci line. Conversely, by 
following all non-parallel Fibonacci lines from the reference plane we reach 
any parallel vertex plane. 

We have shown that all vertex points q\\ such that q± belongs to the Kepler 
Xx model appear in an infinite sequence of parallel planes in the order and 
spacing of a Fibonacci line. The analysis of the ID Fibonacci system given 
in subsection 2.1 now applies to the 3D tiling. The stepwise generation in 
E± of the ID Fibonacci system shown in Fig. 2 can now be converted into 
a stepwise generation in E± of parallel planes, enumerated by the integer 
N: Starting at an LL vertex means starting at a dense plane. The value 
of the perpendicular coordinate y± yields the value of a coordinate tjt ©, 
— 1 < rj < 1 from the center of the triacontahedron along the 5fold axis. The 
explicit relation is 

V (N) = ^y±(N). (10) 

T + Z 

In Fig. 6 we plot the projection of the triacontahedron together with the 
values of y±(N) connected by lines. Each value determines the height of 
a corresponding horizontal section of the triacontahedron. Since the point 
density in the corresponding plane of the tiling is proportional to the area 
of this section, Fig. 6 provides insight into the variation of this density from 
step to step. 

Note that within the Fibonacci window we do not reach the highest absolute 
values of rj. For the additional points in Fig. 6 see section 2.6. 
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Fig. 6 Vertical values of y±(N) ~ 7](N) eq.(10) determine horizontal sections 
in E± of the triacontahedron as windows for planes. The values are connected 
by lines and follow the Fibonacci coding. The total width of the vertical 
window is f ©. The triacontahedron has the vertical diameter 2r ©. 
Numbers in the second row assign planes corresponding to the terraces found 



The bars 1 ,9 ,3 + ,6 + ,ll + mark values of rj for additional low- 



density planes of vertex points. 



2.6 From planes to terraces at the surface. 

We have found in section 2.4 from the bulk model sequences of planes with a 



spacing that agrees with the terrace spacing found in flnj. We now interpret 
terraces at the surface as particular planar terminations from the bulk tiling 
model. From Fig. 2 we have already identified a string in correspondence 
to the observed string of terraces. With the numbers 1 ... 11 in the second 
horizontal row of Fig. 6 we now assign planes and values 771 to eleven planes 
which correspond to the spacing of eleven high or low terraces found in fl~5|| . 
The numbers follow the terraces in a direction into the bulk material. 
The numbers rji given in Table 1 are not unique, but the appearance of the 
finite string puts narrow bounds on their range: Maximal shifts upwards by 
Ar] = ^j(2r — 3) from N = 9 or downwards by A77 = frrjC— 3r + 5) from 
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N = 14 are compatible with the appearance of the string, but of course give 
different values of the density and of the Patterson function. 
As in the patch analysis given in [12|], there appear additional vertex planes in 
E\\ with a spacing scaled by r -1 which do not yield terraces in the experiment. 
This narrow spacing cannot be coded by a single Fibonacci line. With respect 
to the triacontahedron in Fig. 5, it requires a vertical shift Ar] = in E±. 
This shift can be produced by the sum of two vectors of length r(2) pointing 
along the axes 2 and 2' in Fig. 2 respectively. The summed vector (not parallel 
to the 5fold axis) connects points in the triacontahedron only if the initial 
point obeys -^p^ < \q\ < In the selection of points of Fig. 6 the values 
r\ of the corresponding 5 final points are denoted by 1~, 3 + , 6 + , 9~, 11 + . The 
± sign codes an additional vertex plane shifted in E\\ by in units r 
(5) above or below the plane with the fixed number. The positions of these 
planes agree with the patch analysis of fT2]| . All these additional planes have 
a very low density of vertex points. 

In section 3.2 we compute the exact model density of vertex points in the 
planes which decreases with the absolute value of rji. From Fig. 6 it can be 
seen that the planes 3,6,9 in the string have the highest values of \r)i\ and 
hence the lowest density. In the experiment |15j these planes correspond to 
terraces of minimum measured planar size. 



3 Clusters, pentagonal faces and cuts, and cor- 
relations in planes. 

We proceed to an analysis of the more detailed model structure within the 
planes. So far we looked at planes occupied by vertex points from the tiling. 
The full set of atomic positions [|4], [0]] comprises more points in various 
classes. The repetition pattern and variation of the density for parallel planes 
found in section 2 is a general property of the tiling and applies to any set 
of atomic positions within planes perpendicular to a 5fold axis. If atomic 
positions of two different types appear within the same initial fixed plane, 
their repetition pattern and variation of density follows the same pattern as 
found for vertex points, but may propagate from different initial conditions. 
The relative density and the correlation of different types of atomic positions 
sitting within the same plane will then show systematic variations along a 
sequence of planes. 

We consider in this section additional atomic positions from the Bergman 
clusters. The vertex positions of Bergman clusters on the tiling will pro- 
duce points in parallel planes. This holds true in particular for top faces of 
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Bergman clusters which run perpendicular to the chosen 5fold axis. They 
are of particular interest as candidates for the pentagons found in [15|]. Upon 
comparing the height of 6.60 A for the Bergman clusters with the lowest 
terrace spacing of 4.08 A we already conclude that these clusters are cut at 
quasicrystal surfaces. A second larger type of pentagons arises from a top cut 
at the heigth of 4.08 A through five vertices of the Bergman dodecahedra. 
In the full model [|, |llj of AlPdMn these top face pentagons have central 
atoms in a lowered central position. Therefore these pentagon faces would 
produce at their centers the observed holes in the terraces [|15|]. We shall 
examine the correlation of both types of pentagons with vertex points of the 
tiling. 



3.1 Three models for correlations of pentagons and ver- 
tex points. 

We analyze three significant types of planes perpendicular to a 5fold axis 
taken from the bulk model and explore the correlation within these planes. 

(i) In a first model analysis the vertex points of the tiling dominate the 
sequence of perpendicular planes, characterized by the values rji as before. 
Consider the relation of the Bergman top pentagons with respect to these 
planes. As mentioned in the introduction, the centers of the Bergman do- 
decahedra take the positions of those (odd) points of the primitive P-lattice 
which are dropped when going from P to (2F). This has the consequence that 
the Bergman centers can be grouped into shifted planes perpendicular to a 
5fold axis. The planar densities of Bergman clusters in these planes must 
follow the same rules, minima, and maxima as given for vertex positions in 
sections 2.4 and 3.2. 

The height of the Bergman dodecahedra is such that a 5fold vector from the 
plane to the center, passing through a pentagonal face, has the length and 
direction of a typical vector esy, compare Fig. 7. With the chosen enumeration 
of terraces we follow them in the direction of the vector eiy , in the direction 
of the 5fold axis in Fig. 4. Now we look for Bergman clusters with a center 
displaced below and a top pentagon within a fixed plane of vertex points. 
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Fig. 7 A Bergman dodecahedron touches with a top pentagon from below a 
plane of vertex points marked by the top horizontal line. The vector eiy points 
downwards. A typical vector downwards from the plane to the Bergman 
center is e*,\\. A vector like {e\ — e$)\\ — e$\\ runs from a planar pentagonal top 
cut through the dodecahedron, marked by the lower horizontal line, to the 
Bergman center. 

Transforming the vector e§\\ from E\\ to E±, compare section 2.2, one finds 
the coding points for all centers of the Bergman dodecahedra displaced by 
e5_i_ against the direction of the axis 5 of Fig. 5, hence downwards in Fig. 5, 
from the plane coding the vertex points. So the planes with a corresponding 
downwards shift Arj = -^-^ form another Fibonacci sequence that codes the 
Bergman centers. The shifted values 7] 2 for this sequence are given in the 
second column of Table 1. In Fig. 8 we show the two values 771, 772 as functions 
of N — 0, . . . , 24. If the vertical shift from the crosses to the circles equals 
At] = the Bergman top pentagons touch the vertex plane from below. 
This holds true except for the planes N — 6, 11, 14, 19, see Fig. 8. In the 
Fibonacci sequence of planes of vertex points this correspond to an LS vertex, 
but in the planes of Bergman faces to an SL vertex. Therefore no Bergman 
top faces can appear in these planes of vertex points. In the string of planes 
9 ... 19 — > 1 ... 11 this would occur at the planes 11, 14, 19 — > 3, 6, 11. These 
planes carry instead the larger Bergman top cut pentagons, see model (iii). 
Figs. 8, 9 demonstrate that planes formed from atomic positions of different 
type (vertex points or pentagon centers) follow the same Fibonacci propa- 
gation law. Due to systematic shifts in the parameter 77 for different types, 
the corresponding densities, which are functions of this parameter, propagate 
differently even for atomic positions of different type within the same initial 
plane. 
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012 345 67 89 10 11 12 13 14 15 16 17 18 19 20 2122 23 24 



Fig. 8. The values y±(N) ~ T](N) for N = 0, . . . , 24 in model (i) determine 
horizontal sections of the triacontahedron as windows for planes of vertex 
points (771 , crosses) or of centers for Bergman top pentagons (772, circles). 




012 34 56 789 1011 12 13 14 15 16 17 18 19 20 21 22 23 24 



Fig. 9. The values y±(N) ~ r)(N) for N = 0, . . . , 24 in model (ii) determine 
horizontal sections of the triacontahedron as windows for planes of centers 
for Bergman top pentagons (771, crosses) or of vertex points (r/ 3 , circles). 
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(ii) In the second model analysis of the planes determined by rjx we assume 
that they are dominated by top pentagons of Bergman dodecahedra located 
below the plane, compare Fig. 9. The additional presence of vertex points on 
these planes is now coded by a positive shift Ar/ = In Table 1 we give 
the new values 773 for these vertex planes and show them in Fig. 9. Additional 
vertex points appear except in the selected planes 9, 17 — > 1, 9. 

(iii) In the third model analysis we assume that the planes determined by 
r/i are dominated by pentagons, scaled by r and corresponding to top cuts 
through Bergman dodecahedra below the plane. A typical vector from the 
plane for this cut to the center of the Bergman cluster in the notation of 
section 2.2 is (— e§\\ + eiy) — e$\\ indicated in Fig. 7. The edge size (2) = 4.78 
A of the pentagons would be in line with the holes observed in UTS] , and 



it would also lead to a central hole. Again we ask about the presence of 
additional vertex points. From the relative position of these pentagons to 
the vertex planes as shown in Fig. 7 we deduce in E± a shift A77 = ^q^. At 
N = this shift modulo the window size becomes which then generates 
the values r/ 3 given in Table 1 and shown in Fig. 9. Vertex points can occur 
only if the relative shift from crosses to circles is Arj = r] 3 — rji = ^r^- This 
happens only in the selected planes 9, 17 — > 1, 9, otherwise there are no vertex 
points within these planes. A further shift analysis shows that any plane 
containing Bergman top pentagons cannot contain Bergman cut pentagons 
and vice versa. In model (iii) there appear Bergman top cut pentagons in 
densest planes. A closer inspection shows that these pentagons share vertices 
and form an almost connected graph. 

The three model cases considered yield three interpretations of bulk planes 



as terminations for the terraces observed in flyj, all in line with the relative 
spacing. In cases (i,ii) the planes carry vertex points or pentagons corre- 
sponding to the faces of Bergmann dodecahedra. In case (iii) the planes 
carry the larger cut pentagons and almost no vertex points. 

3.2 Planar density of atomic positions. 

We compute the exact area F(r}) of a planar section of the triacontahedron 
as a function of 77, — 1 < 77 < 1. This function is proportional to the exact 
density of vertex points D(rj) in the plane coded by this section. The result 
is 

0< \V\ (11) 
t + 2 

F( V ) = (r + 2)- 3 / 2 [lOr] , 
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The function F(rj) is plotted in Fig. 10. The maximum is -F(O) = 10r(r + 
2)-3/2 _ 2.3511, the maximum value of 77 for the Fibonacci sequence of planes 
is 

|77| = r 2 (r + 2)- 1 = 0.7236. (12) 

These values are marked by vertical and horizontal lines in Fig. 10. 
The density can be converted into the absolute density of vertex points by 
considering the densest planes with a triangle pattern and its vertices: In 
the triangle pattern, each triangle contributes, because of the sum \ 2ir of its 
angles, a weight \ to the number of vertex points. In terms of the short edge 
length s, the area fi, fi and relative frequency i>i, u 2 of the large and small 
triangle are 

h = s 2T - v^+2, j 2 = r" 1 /!, v 2 = r- 1 ^, (13) 

These expressions yield for the density of vertex points (number of vertices 
per unit area) the value 



D(0) 



2r 2 



2/ir + 2 



(r + 2) 



3/2 



In the present model we put for the short edge length s = r (2) = r 
(5) and (5)= 4.56 A eq.(l) to obtain the model value 

D max , {i) = D(0) = 12.6 • 10' 3 A" 2 . 
The density of vertex points in a plane for fixed 77 is now computed as 

F(rj) 



D(v) = D(0) 



F(oy 



(14) 



/r+2 



(15) 



(16) 
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The lowest density in a Fibonacci sequence of planes from eqs. (14,15) is 

D mm = ^-D(0) = 0.309(LD(0). (17) 
It 

These values of the density refer to vertex positions in the tiling. All other 
atomic positions on the tiles will of course lead to other specific atomic den- 
sities and correlations. 

The maximum density of vertex points, associated with model (i), is given 
in eq. (16). The same maximal density applies to the centers of Bergman 
top faces in model (ii). Each pentagon contributes 5 vertex positions which 
would yield a factor 5 for the pentagon vertex density. For model (iii) with 
the larger pentagons, the maximum density of centers is still the same. For 
the density of pentagon vertices one should not multiply by a factor 5: It 
turns out that these larger pentagons in a dense plane can share vertices. An 
exact computation of the maximum density for the vertices of large pentagons 
yields (proof omitted) 

It + 4 

D ma x,(iii) Dmax,i 3.6180Z) maa;) (j) . (18) 

From the density we can compute the average distance by comparison for 
example with a tiling by equilateral triangles. For such a tiling, the density 
of points is related to the edge length t by 

Dsit) = r2 73 (19) 

If we equate this density with the expression eq.(15) found in the dense vertex 
planes, we obtain an equivalent distance t eq = 9.5 A . 

In columns 5-7 of Table 1 we give F for the three values of 77. Note that, in 
view of the three models discussed in subsection 3.1 , the three values of F 
in a row do not always refer to the same plane. 



21 



0.5 




V 

Fig. 10 The function ^(M) for a plane fixed by 77 > is in E± the area of a 
planar section of the triacontahedron. The horizontal line marks the lowest 
value of F, the vertical line the highest value of \rj\ in a Fibonacci sequence 
of planes. In E\\, F(\r)\) is proportional to the density of points or pentagon 
centers within a plane. 
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Table 1: The values 771,772,7/3 as functions of N = 0, . . . , 24 code inside 
the triacontahedron Fibonacci sequences of planes perpendicular to a 5fold 
axis. These values are associated alternatively with vertex points, pentagonal 
faces and cuts of Bergman dodecahedra. Columns 5 — 7 give the value of 
F(r]i), equal to the area of planar sections through the triacontahedron and 
proportional to the relative density of points. The rows 9 ... 19 are put in 
correspondence with the terraces 1 ... 11 of flyj| . 



N 


'/l 


Tin 

'12 


Tin 

'13 


vn) 


F(rin) 

\'I2) 


vis) 


n 

W 


-0 1 708 


-f) 4479 


1056 


2 351 1 


1 9021 


2 351 1 


1 

J. 


7236 


4472 


-f) 4472 


7265 

yj . § — u • j 


1 9021 


1 9021 


2 


1 708 


-0 1 056 


4472 


2 351 1 


2 351 1 


1 9021 


Q 

O 






U.1U-JU 


9 DRQ1 


1 0^41 


9 3^1 1 


A 

1 


51 25 


2361 

W . — ■ )\) ± 


-0 6584 


1 6746 


2 3261 


1 0541 


•■J 


-0 0403 


-0 31 67 


2361 


2 351 1 


2 2260 


9 3961 


6 


-0 5931 


5777 


-0 31 67 


1 3507 


1 4169 


9 9960 


7 


301 3 


0249 


5777 


2 251 


2 351 1 


1 41 69 


8 


-0.2515 


-0.5279 


0.0249 


2.3129 


1.6164 


2.3511 


9 


0.6430 


0.3666 


-0.5279 


1.1269 


2.1259 


1.6164 


10 


0.0902 


-0.1862 


0.3666 


2.3511 


2.3497 


2.1259 


11 


-0.4626 


0.7082 


-0.1862 


1.8512 


0.8067 


2.3497 


12 


0.4318 


0.1554 


0.7082 


1.9508 


2.3511 


0.8067 


13 


-0.1210 


-0.3974 


0.1554 


2.3511 


2.0495 


2.3511 


14 


-0.6738 


0.4971 


-0.3974 


0.9796 


1.7311 


2.0495 


15 


0.2207 


-0.0557 


0.4971 


2.3365 


2.3511 


1.7311 


16 


-0.3321 


-0.6085 


-0.0557 


2.1982 


1.2835 


2.3511 


17 


0.5623 


0.2859 


-0.6085 


1.4800 


2.2733 


1.2835 


18 


0.0095 


-0.2669 


0.2859 


2.3511 


2.2969 


2.2733 


19 


-0.5433 


0.6276 


-0.2669 


1.5565 


1.1980 


2.2969 


20 


0.3512 


0.0748 


0.6276 


2.1600 


2.3511 


1.1980 


21 


-0.2016 


-0.4780 


0.0748 


2.3456 


1.7986 


2.3511 


22 


0.6928 


0.4164 


-0.4780 


0.8851 


1.9966 


1.7986 


23 


0.1400 


-0.1364 


0.4164 


2.3511 


2.3511 


1.9966 


24 


-0.4128 


-0.6892 


-0.1364 


2.0070 


0.9033 


2.3511 



The three models (i,ii,iii) given in subsection 3.1 yield different density values 
F. In model (i), the density of vertex points and centers of Bergman top faces 
is given by F(r)i) and, with exceptions, by F(ri 2 ). In the plane 16 — > 8, the 
density of Bergman faces becomes D = 6.8 10~ 3 A -2 , respectively. In model 
(ii), the density of Bergman top faces and vertex points is given by ^(771) 
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and with exceptions by F(r] 3 ). respectively. In model (iii), the density of 
Bergman top cuts is again given by F(r)i). 

The density of pentagonal holes in the experimental data [T5|] has approxi- 
mately |12| the value D = 4.2 1CT 3 A~ 2 . This value would favour model (i) 
with Bergman top faces in planes dominated by vertex points. 

3.3 Patterson analysis in planes. 

Consider the general Patterson function P(x\\) of a quasiperiodic point set 
with a window W at each point of the lattice A. Let v\\ be a shift vector 
projected from a lattice vector v. Let x( x ±) be the characteristic function of 
the window W for the lattice points projected to E±. 
Prop 6: The Patterson function in E\\ at x\\ — x — x± is given by 



P tota '(x,|) = <*(*|| ~v\\) X(x±)x(x ± - v ± )dx ± = J2 S(x\\ - v\\)P(v ± ). 

v<=A JW v&A 

(20) 



We shall put P(v\\) = P(v±). A Patterson analysis within a plane perpen- 
dicular to the 5fold axis 5 with a fixed shift vector v\\ parallel to this plane 
involves the following notions in E± : Consider the corresponding planar sec- 
tion of the triacontahedron with an intersection at r] r (5) along the 5fold axis 
5. The values of r] for the 11 selected planes are given in Table 1 both for 
the vertex points and for the centers of Bergman dodecahedra touching the 
vertex plane from below. The area of the planar section is proportional to 
the density D{rj) of vertex points. The terraces 2,5,10 in this interpretation 
are dense vertex planes. The terraces 2 and 7 have the highest, the terrace 
6 has the lowest density of Bergman faces. 

The Patterson function from Prop. 4 is computed as follows: Shift the planar 
section by the vector v± parallel to itself. Compute the area of the intersec- 
tion of the shifted and unshifted section. This area is proportional to the 
value of the Patterson function in E\\ at the point v\\ . The value of P{0) is 
proportional to the density of points. We may also normalize by plotting 



In a circle approximation we proceed as follows: 

We use the exact area F(rj) of a planar section of the triacontahedron as a 
function of r). 

For any fixed value of r/ we compute the radius r{vj) of a circle with the same 
the section, hence 



P(v\\)/P(0). 




(21) 
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This amounts to replacing the triacontahedron by a rotational surface whose 
circular areas (and hence density values) for any 77 are equal to the ones 
of the triacontahedron. Then we approximate the Patterson function by 
using these circles instead of the polygonal sections as functions of 77. In this 
approximation, the Patterson function only depends on 
The Patterson function P has a smooth behaviour in terms of the two vari- 
ables 7] and v±. This expression is shown in Fig. 11. For v± = it reduces to 
the function F(rj) and for v± = 2r(0) it must go to zero. 




Fig. 11 The Patterson function P(r),v±). For v± = it reduces to F{rj). 

For the plane number 16 — > 8, we give in Figs. 12-14 the values of the Pat- 
terson function P for 771,772,%, represented by areas of circles, as a function 
of the 11 points v\\ selected in a plane as in fllBl . 
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Figs. 12-14. The Patterson function P corresponding to the values N = 16 — > 
8 for the values 771,772,^3 from Table 1, represented by the area of circles in 
a plane. The upper right-hand circle stands for : v\\ — 0, the first point I' 
to its left is at a distance \v\\\ = 7.78 A the others are labelled by roman 
numbers 7 ... X in correspondence to JTofl. 
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In Figs. 15-17 we present the same values as a function of the 10 roman 
numbers which label peak positions v\\. It can be seen that P(r} 2 ) yields the 
lowest density and the strongest relative variation, P(rji) and P(fj3) are very 
similar, but the latter yields the highest density. 




I' II IV VI VIII X 




Figs. 15-17. Patterson function P as in the previous Figure as a function 
of the roman indices P,I,...X for the three values f]i,f]2,V3 from top to 
bottom. 
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We can now compare the three different models (i,ii,iii) of the terrace struc- 
ture given in subsection 3.1 in terms of the Patterson data for terrace 8. 
Model (i) yields the lowest value of the Patterson data for the Bergman top 
faces, along with a larger value for the vertex points. Models (ii) and (iii) 
give the same larger values of the Patterson function for the small or the 
large top pentagons, but differ in the presence or absence of vertex points 
respectively. All three models are compatible with the qualitative set of 



experimental Patterson data [15 



4 Conclusion. 

We analyze a tiling model for the surface structure of i-AlQ 8 Pd 2 3Mn 9 qua- 
sicrystals perpendicular to a 5fold axis. The surfaces are interpreted as ter- 
minations in atomic planes from the bulk tiling model decorated by Bergman 
clusters. The model is analyzed in a window approach. The quasiperiodic 
relations between planes occupied by atoms are rigorously and quantitatively 
derived from the bulk model. Sequences of planes are shown to be generated 
and connected by infinite Fibonacci lines running through the tiling. Both 
the density and the Patterson data for different types of atomic positions are 
computed in closed form and shown to vary strongly between the planes. 

In the planes we consider vertex positions of the tiling and pentagonal posi- 
tions arising from faces or top cuts through the Bergman clusters. Changes 
of their correlation between the planes are analyzed in three alternative mod- 
els. The pentagonal holes observed in [|15|] admit an interpretation in terms 
of these pentagons. If they arise from the Bergman top pentagons, their 
observed larger size must reflect a local reconstruction of the surface. If they 
arise from the Bergman top cut pentagons, they agree with the observed size. 
The geometry, the spacing, the density and the Patterson functions are com- 
puted from the bulk model. The sequence and spacing of terraces and the 
available Patterson data from |15| are well reproduced. The observed Fi- 
bonacci string of planes yields information on the density of atomic posi- 
tions. The observed size of the terraces shows some correlation with the 
model structure in the planes. 

More detailed experimental studies of the terrace structure in i-quasicrystals 
are suggested. 
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